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Abstract: The time diffusive boundary layer from rest for Taylor series uniform outer flow is constructed by 
integration of the complementary error function of the diffusive similarity variable . The same erfc solution is 
also easily “Stokes ” transformed to solve harmonic outer flow from rest. 

If the outer flow vector varies weakly along the boundary, the perturbation pressure gradient accelerates 
the slowest fluid nearest the wall the most to disproportionately vary the shear stress. So centripetal 
acceleration causes secondary crossflow inside the boundary layer with strong wall shear towards the center of 
curx’ature. But longitudinal acceleration also implies an inflow towards the wall which thins the outer boundary 
layer with a weak further increase of the wall shear, but the strongest perturbation steady streaming. 

Simple new particular solutions for these two perturbations are easily constructed in terms of products of 
integrals and derivatives of any primary diffusive solution. For an outer flow as a series in the square root of 
time, all homogeneous time coefficients remain just iterated error functions. Each systolic pulse in the aortic 
arch was considered as a Taylor series flow from rest to calculate the wall shear vectors. 

When the outer flow oscillates forever more, its primary diffusive boundary layer asymptotes to the Stokes 
oscillatory exponential decay with distance from the wall. The particular perturbations are exactly evaluated 
and also confined near the wall but with mean slip. The mean slip homogenous perturbations diffuse outside the 
Stokes layer into steady streaming as complementary error functions with inverse time correction functions. 
Extended Taylor series computations provide more detail of the perturbation transients. 
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I. Introduction 

A theoretical study [1] was undertaken of the 
boundary layer in the aortic arch with a view to 
possibly explaining the patchy occurrence of 
atherosclerotic lesions on the arch walls, for 
instance due to variation blood particle motion near 
the wall [2], Experimental study of the boundary 
layer is complicated by the flexibility of the aortic 
wall and the difficulty of in vivo examination, 
requiring the sacrifice of large animals such as 
horses. 

The local diffusive boundary layer was perturbed 
for pressure gradients and inflow due to the helical 
curvature of the aortic arch and the stream wise 
variation of that curvature’ strength and plane. The 
maximum bend is tight and it was not necessary to 
resort to the inappropriate small curvature 
approximation commonly used. In fact the 
numerical results were quite different than the low 
curvature ones, with secondary flow on the inside of 
the bend much intensified by the stronger potential 
vortex core flow there. This indicated the secondary 
flow would focus into a jet normal to the wall at the 
very inside of the bend. The jet separates the 
boundary layer and injects vorticity into the flow as 


a whole. This intrusion of the boundary layer into 
the bulk flow is seen in the pulsed Doppler 
ultrasound profiles of Peronneau [3]. Downstream 
in a light coil, this core vorticity would eventually 
transform the flow into the high Dean number 
steady flow of greater axial velocity near the outside 
of the bend. 

Some correlations were found between wall shear 
variations and the patterns of disease which is 
actually different from that of dye uptake and other 
experimental patterns. Unfortunately no absolutely 
compelling matches were found and the entire 
physiological localisation patterns and possible 
mechanisms need closer scrutiny rather than further 
flow studies. 

During this work simple ways were found to 
construct homogeneous solutions for the diffusion 
(heat) equation and first order solutions for fluid 
spatial acceleration perturbations of it. These might 
also be adapted to solution and perturbation of other 
linear pde’s. Here they will be presented for the 
generalised fluid dynamic problem of perturbing 
starting and oscillatory boundary layers for spatial 
acceleration of the coreflow. 
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Figure 1 Illustrative Geometry 


II. Perturbed Diffusive Boundary layer Equations 

As a simple illustrative case consider the inviscid incompressible streamline velocity G(z)W'(X) parallel to 
and small distance Y below the centerline of the top wall of a slightly tapered square duct, as in Figure 1. Here 
G(t) is a non-dimensional 0(1) pulsatory/oscillatory time factor starting at time t= 0. Then the axial X gradient 
of pressure p due to temporal and advective accelerations resp. is the inviscid 

d xP = -p (Wd r G + G 2 Wd x W ) (1) 

with p the density. With v the kinematic viscosity, the 2D boundary layer equation for the centerline velocity w 
at small distance Y below the top wall is 


d T w + wd x w + udyW = -d K p/p + vd 2 y w (2) where d>« + d v iv=0 (3) 
w>=0 at Y=t= 0 and w tends to GW as Y tends to oo. 

Let t= t IT , T a characteristic time such as the angular period of a harmonic G, and let y=17 V v/\ the 
characteristic diffusion distance scale in time T 

Then d r - vd 2 y becomes D/T where D= d, - d 2 y is the Diffusion operator 

When the VFd fi dominates G 2 Wd x W, the approximate primary diffusive O’th order solution near r=0 is velocity 
Wg 0 (y,t) with non-dimensional deficit or interference g= G- go such that 


d,g-d 2 ,g=D[g]=0 g=G(t) aty=0 (4) 

Now the pressure gradient due to the acceleration of the coreflow parallel to the wall d x p=-p G 2 Wd x W 
will perturb this boundary layer more than the correction for the acceleration in the boundary layer p wd x w so 
the net “planar/pressure” l’st order correction is TTFdAWvi where Td x W«\ 


D[ Vl ]= G 1 -gf=2Gg-g 1 (5) 

The first term dominates and the net peaks at G 2 at the wall which of course restrains the response. Thus very 
strong shear is to be expected from it as of interest in atherogenesis. For instance, if the duct is also curved with 
centerline radius R as in Fig 1, the inviscid flow has a lateral centripetal pressure gradient dzP=-pG 2 W 2 /R so 
that in the boundary layer a secondary lateral flow of v= T W 2 Vi /R starts towards the center of curvature. 

For a longitudinal acceleration the walls are converging towards the center axis flow but the normal component 
u must decline to nil at the wall ie from (3) u= -dxWJo 1 go dF which increases at first as the square of the 
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distance from the wall where g„=0 until linearly outside the boundary layer. So there is also a positive l’st order 
correction to g 0 for - udyW “advection in” or “indraft” as TWdxW w i 

where D\w,\= djgJo’ g 0 dT= d y g o yg„dy (6) 

For a typical (starting) boundary layer profile such as g 0 = Gerf r|, drg,, JoVo dr| peaks at about G 2 /10 at r) ~1 
(where 2 Gg-g 2 ~.3) before vanishing at infinity outside the boundary layer. So the W\ response to this forcing 
can be expected to be about .1 Vi‘s and this is indeed borne out for the shear perturbations in both starting and 
oscillating coreflows. But anticipating section IV, comparing r| moments about the boundary for the same G, 
the factor is about .3. Now 

Wg 0 d.V = W(G-g)dy= {yG-jVgdy } so I)\w x \= d v g {-yG+Vgdy} (7) 

where the first term is the inviscid inflow dominating away from the wall, though near the wall it is cancelled to 
leave the nety 2 increase. So then the first order longitudinal acceleration correction terms to w are TWdxW 
(Vi+w'i), small vs W as assumed if TdxW« 1. 

Thus there are two basic and very different perturbed diffusion equations to be solved in all first order 
perturbations of a time diffusive boundary layer for spatial acceleration effects. Previously wq+Vi longitudinal or 
Vi transverse have been laboriously solved for specific G(t) eg impulsive, linear and quadratic or at t= co 
harmonic [4] without interconnection or comparative understanding. Here general (particular) solutions for the 
w i and Vi are presented and then the complete solutions found for any series in 7, and asymptotic solutions for a 
simple harmonic. 


III. Perturbation Solutions for General Time Variation 


Since D is a linear differential operator with homogeneous solution g: D[g]=0 , then d,g or d v g are further 
solutions G[d,g]=0 and £>[d v g]=0. Likewise the indefinite integral in t generates a solution, or the definite 
integral fj gdt if g=0 at t=c, and //gdy if g=d y g=0 at y=c. To generalize the inversion e ry f c v e' ly gdy of the 
basic differential factor d y +r commutes with D provided e' 1 g & e' ly d y g vanish at y=c. (This provides one way 
of generating the g for harmonic G(t) starting at f=0 in section IV). 

To find particular solutions of the above perturbed diffusion problems, simply note that D\gE(t)\=g<\,H and that 
if h is another homogeneous D[h]= 0 solution, then />[g/i |=-2 d,g d y h These immediately provide a particular 
solution to the simpler pressure problem Vi= Vi p - Vih the difference of particular and homogeneous parts such 
that 


D\v lh \=() /7[V| P J=2Gg-g 2 v lp =2g/ 0 t Gd/ + Vi[jJgdyf (8) 
Then aty=0, Vi p =Vi h =2G/ 0 t Gdf+ l /2[/„ o gdy] 2 (9) and d y Vi p =2d y g/ 0 t Gdf -G/ 0 ”gdy (10) 

For ivi (7): £>[wi p ]=d v g {-yG+J 0 y gdy} 

Trying -y d y g I a Gdt gives the dominant first term plus 2d 2 v g/ 0 f Gdt or 2d,g/ 0 f Gdt which like the last term d y g 
/o y gdy} = d v g {/o'gdy+/./gdy} is easily solved by the above properties to get 

ivi p = -y d v g/ 0 f Gdr -2d, gfj dt Ij Gdt +d y g ffgdy -Vig fj fj gdy (11) 

aty=0 w lp =w lh = -2d, Gfj dtfjGdt +d v g/ 0 7ogdy - ViG/f/Jgdy £>[w lh ]=0 (12) 

d v >v lp = -dygfa'Gdl -2d 2 , v g/ 0 t dtfjGdt +d, G fjffgdy - Vi d v g fjfj gdy- 1/2 G/f gdy (13) 

there are no common elements to Vi p and iV| P (because their fluid dynamic origins are distinct) so they are the 
simplest basis functions to solve versus any linear combination such as Vi p + iv, p 

IV. Diffusive and Perturbation Solutions for Root Powers of Time 

From the fundamental solution (?o=erfc r| of the diffusion equation D[Qq\= 0, where the diffusive similarity 
variable r\=Viyl jt , the same d y and f~ y dy operators can construct series of solutions Q„ for g=G(t ) =7" at 
y=0 as in Table 1. 
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TABLE I D[F„]=0 solved by Q n and P n of form F n =^t"f n (r\) F n =d, F ll+] r\=Viy/jt 
With q and p both solving f(r\) such that f n -Af ’ n+ i nf n -x\f n ‘-Vif n ”=0 
so Vmf„= r|. 2 Define y n =T[w+ Vi ] 


;z 

jt" 

q„( ti) 

?«(0)=(-l)7y„ 

Pn( 0) 

Pn(f])= Vi lim q n (r[) riY- 

4 

t 2 

16i 4 erfc(ri) 

V 2 

V 2 

'A +2ri 2 + 2rj 4 /3 

3 

jt 3 

-8i 3 erfc(rp 

-4/3^71 

0 

2t!+4ti73 

2 

t 

4i 2 erfc(ri) 

1 

1 

2t[+\ 

1 


-2iierfc(ri) 

-2/Vji 

0 

2ti 

0 

1 

erfc(ri) 

1 

1 

1 

-1 

i/Vt 

-e^ 2, jn 

-i/Vjt 

0 

e z2 dz =p.i 

-2 

r 1 

T^e^ 2 Nn 

0 

1 

l-il P- 1 

-3 

t -m 

('A-r\ 2 )e~ r]2 /^k 

'A Hn 

0 

(r[-'/i)p.\ -r| 

-4 

f 2 

(T| 3 —3rj/2 )e-" 2 Nn 

0 

-1 

(-Tl 3 +3i 1 /2) i ;. 1 + (ri 2 -l) 


Notes: The multiplier of -e 1,2 /Vtc is the same as the multiplier of p. t = 2DawscnF(r|) The negative index 
P. n = d""\ (2Dawson F(r|)/) so p. n = (V 2)"' 2 d 11 "' 1 DawsonF(ri) which alternates in non-zero value at j’=?p() 
with the q.„ to allow solution for G(t) as asymptotic series in inverse powers of large time. The p n series 
breaks between f jp.idy divergent and d v /7 0 =0, and for positive index the p„ diverge at infinity. 


These solutions of the general form F„=jt"f„( r|) so d, =d ri / 2Vt dy=2vt dr). 

D= d t -Va\ t % - 'At -'8\ so D[Fj\ =V 2 ^ n - 2 (nf,-r\f n ‘-V 2 f„”) =0 

a second order ode with two independent set of solutions being the Q„=st"q n (, T|) and P n =yt"p n (r\) 

tabulated above. Since all/„(ri) are series of r| derivatives, this o.d.e gives the recursion relation. 

V 2 nf„= r)/„.i +/„.2 useful for computing. 

For non-negative n, q„ =(-2)" i„erfc(ri)=(-2)" i„ the iterated complementary error functions as in Figure 2 
where d n i„= -i„.i. All q„ or i„ vanish at infinite r|. 

So for G(t) = E a n jt n (with fjGdt= E 2 a n jt n+2 / /i+2) g= E a n y„ (2V/)"i« where the y„= T(Vm+\) 
normalize the 2 n i„ as in Figure 2. Note y„/ y „.2 = 'An. 

Remember i„ + i=Jr, K d n i„ so the primary shear is d,g= - Ea n y„ (2\f)"" 1 i„. ] . Thus with a pulsatile flow for 
example a parabolic pulse, the shear reverses before the coreflow. Then from (8) and (9) 

Vi=EE «„a s (2 + ' +2 [y„i„/2'(.s'+2) +1/ 2 y„i, i+ iy s i s+1 -y„ +s+2 i, i+s+2 {4/(s+2)+ Vi y„y J y n+ iy s+ i}] (14) 

As anticipated the Vi p or first two terms in the bracket are positive with negative shear reflecting the increase 
of the forcing towards the boundary. So the positive shear of the negative {[homogeneous solution must 
prevail, since the overall shear should be positive. Physically the strong slip of the particular solution at the 
boundary from the forcing generates strong homogeneous shear, dominating the weak reverse shear from the 
decrease of the particular solution away from the boundary. So from (10) and (14) 

aty=0 d v V|=EE« n « s (\t)' ,+1+ ’[-4y„/ y„.,(.v+2) - y„/y„ + i+y„ +s+2 /y„ +s+ i{4/(s+2)+ 1 /2y, ! y s / y„ + iy s+ i}] (15) 
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Comparing the two particular wall shears for n=s, 4y„ +1 / y n .\ =2(h+1) vs. (h+2) so the first particular term 
dominates with increasing n from equality at n= 0. Regarding Figure 2 as relative velocity profiles it is clear that 
the shear of the homongeneous normalised i 2n +2 dominates the first particular and so the sum even though it is 
less than the shear of the square of i„+i in the second terms. 

Evaluating the w i perturbation requires further 

If gdy=La n y n {2jt) n+ \ +l d n g=-S« n y„( 2 Vf)"i„. 1 , d,G=Z 1/2 najt "- 1 , 

/o7o* Gdt=Z4aft n+4 / (b+2)(b+4), f 0 l g=UfJgdy=Z a n y n (2jt) n+ \ +2 , d lg =d 2 y g=Za n y n (2ltr \. 2 , 
d\g = d 3 y g=-I.a n y„ (22t) n \. x fc ff gdy=Z2a n y n ftf +3 /y n+l / h+3 , 

I-JIJ gdy=Za n y„ {2jt) n+ \ +1 jJ fj gdy=Y2a n (2 Vf)" +2 / h +2 
so that 


Wif&Z a„a s (V0 n+s+2 [Ti2' !+1 Y n i B .i/(s+2)-2" +1 Y H i n . 2 /(s+2)( S +4)-2Vn-iY s / y s+ iG+3)-2" +I+1 Y s i s Y„i„ +2 1 (16) 

The first term dominant at average r| vanishes at the boundary leaving the negative corrections terms with 
positive shear. But the positive shear of the dominant term means the negative shear of the positive 
homogeneous terms is overcome, w 1 = w i p - w i h Here the wall shear is dominated by the particular solution 
increasing with distance from the boundary like the forcing just overcoming the weak homogeneous shear from 
the weak slip of the particular solution . the opposite situation from Vj. 

M , ip= w lh = EE a n a s VC’ +s+2 [ -4h/(s+2)(s+4) -2y„y 5 /y„.,y v (v+3) -Vzy, /y„ +2 ] aty=0 (17) 

Wih= -EE « n fl s (270" + ' +2 |4H/(s'+2)(v+4)+2y „yj y„.iy s+ ,(v+3)+ 1 / h+2] y„ +s+2 i n+s+2 (18) 

Aty=0 d v H’in= E Er/„a s (2\V)' !+ ' +l [ “ “ ] y„ +I+2 /y„ +s+2 (19) 

dyWip =EE a n aft " +s+1 [2y,/y„.|(v+2)+8y ( /y„.V(v+2)(v+4)+HY s /y (+ ,/(v+3)+y,/y„. 1 /(s+2)+ 1 / 2 y,/y„ + i] (20) 
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This has a very high degree of reduction by the homogenous shear so the net positive shear is as expected about 
1/10 th of Vi’s declining similarly with n. So in the combined axial acceleration problem Vj will totally dominate 
the shear, as anticipated due to the wq forcing peaking low in the middle of the primary boundary layer. 

The aortic systolic flow pulse was represented by a quartic polynomial in t by and the above diffusion 
approach was extended to higher order and to the aorta entry by beginning diffusion at the f>0 when the Oseen 
.35 times the inviscid flow left this leading edge of the wall. The 2 first order perturbation and the linear r terms 
then agreed to +/- 20% with perturbations of the Blasius solution. [1] 

But the sharp curvature in the aortic arch makes the Vi and higher secondary flow perturbations very strong, 
indeed too strong at the inside of the strongest bend where the secondary flows converge to a stagnation point, 
and so thicken the boundary layer that it separates. Actually this begins just downstream of the maximum bend 
where the axial coreflow is also decelerating longitudinally on the inside of the bend. The separation of the 
boundary layer in the computations was seen in ultrasound measurements in humans by Perroneau[3]. 


V. Stokes Transformation, Layer, and Steady Streaming 

The Stokes 'layer' S= e" +y r , d v S=SVr £>[S]=0 solves the diffusion equation with complex Vr 

If any solution at all D[(?]=0 is expressed in terms of t and r| ie as Q(t,x\) and ‘ again denotes differentiation wrt 
the second (similarity) variable r\=Viyl jt 

d t Q- L /zr\Q , /t- 1 AQ'”/t=0 again. Now if £=ri+V rt then d,^=-q/r+V 2 \)/r=-^ /2t+jr/t so 

D[SQ(t,Q]=D[S]Q+SD[Q\-2d y Sd y Q=S{d l Q -ViC, Q’/t-'AQ"lt)+ Q’jr/t }-2rSQ ’d v t]=SQ’ jr/t-2SQ VNr/t =0 

So “Stokes transforming’ ’ any solution Q(t, tj) may give one other solution Q r =S Q(t, rj+\rt). Transforming 
again at -\V reverts to the original Q(t, rj) or at jr just transforms once at 2 Vr . Since ( rj+krt ) 2 = jf+yjr +rt , 

Sexp -(r\+jrt) 2 = exp-r| 2 , the (f of negative index aren’t changed 1 but for all the Q n of positive index 
Q„ r =e" +J T2Vf)"i„erfc(r|+v i rt) are distinct solutions. A less general, more arduous construction is 
Q„, :^ r I Je- r ‘ r Q„ A , r dy 

Note erfc(z) + erfc(—z) =2 so the solution for G(t)= e rt complex exponential starting at t= 0 is 2 
g-MQoAr +Qo,Ar)= 'AV‘W" erfc(q+\W)+ e'^' erfc(r|-\W)}. (21) 

y integrations of Q 0 , Vr=e r /~ y dy e' y '‘ e _T|2 / Vji t for the particular solutions are possible by parts. 

Eg. Jj gdy= 'A(Qn, ,r -0o,-VrWr , JJ JJ g<\y='MQ % , r +Q 0 ,.^)/r+ Q 0 Nr 

So gdy= e" erf (\'W)/ Vr as in Fig 3 but upon squaring this has defeated exact homogeneous solution Vih . 
But asymptotic ones will be considered after clarifying the asymptotics of g and how real g relates to the 
complex. For large | rt | »tj or t » Viy! \ r \ erfc(z)~ e -z2 /z\% (l-!/ 2 /z 2 +..) for large z [5] . Then 

g ^N- Sry -Q,Jr-QJr-+...(22) 

where the second term as t 1 will asymptotically dominate for all r with a negative real part. As Q. 2 =d, Qn = 
d,erfer| it represents the languishing as t _1 at diffusing distance rj of some effective net g from small t. For 
positive real part of r 2 or exponential growth g~S~ e rt " r ' the Stokes layer has an oscillatory component in y 
if G has one in t. For pure oscillatory imaginary r at k'th harmonic r=ik. Choosing and henceforth implying 
the Real part by a cos t (impulsive) start avoids the imaginary Q.i/r so the correction is - Q.Jr 2 as t " 2 which 
doesn’t complicate the study of the asymptote to steady streaming. Then at large t 

g~e"' y '‘ +Q- 4 ~e yH2 cos(t-yN 2 )+Q. 4 , d y g=-e' y ' 2 cos(f-y/V2+%7t)+g. s> J,"gdy=-cos(f- 'Anj+Q.i (23) 

1 though transforming Dawson’s integral solves the unbounded G(t)= i t erf ijrt. 

~ Recursion gives iierf'c(z)= -zerfc(z)+ Q.\ so even in z g= (Qi,. r - Qi,N 4Vr solves G(t)= te" from t= 0 
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As in Fig 3 the transient between J , 0 gdy= Re{-e' ( erfVi/)/Vi} and -cos(t-Vm) can alternatively be uniformly 
approximated by the sum of two exponential decays . 18 e'' 4 ( +. 5 e‘ 4( . The square of this accurate / x ° gdy 
would give Vn, at y=0 in complex exponentials again, which can be solved 



by the above with complex r. Including these decays would give small corrections to the following coefficients 
of Q. 2 or t _1 asymptotes which arise entirely in the homogeneous problems. With the asymptotic g, the 
solution Vi is straightforward to determine 

Vi~2 e' yi ' 1 cos(t-y/j2) sin/ + Vi e' y:1 cos 2 (/-y/\ 2 - 'An) - 5 e' 3 sin( 2 /-y )/8 - 'A erfc77 - 5 Q . 2 /8 (23) 

The square of the y integral of the primary Stokes oscillation in y and t thus creates a mean particular flow at the 
wall whose homogeneous reaction diffuses to infinity as a Q 0 =erfc^. The second harmonic homogeneous terms 
carries a Q_ 2 or /" correction from (21). Because the particular solutions asymptote faster as / " there is no 
point in using the exact ones, but those do show the asymptotic particular solutions below do not need to be 
corrected for initial values. From (10) and (23) 

d y Vi~ , /4\2+(5/4\2-3/2)sin(2/+'/47t) - l / 4 Q.,-%Q.,=. 353 +. 267 sin( 2 /+'/ 47 t) +W V/tt - 5 / 16^1 nt 3 (24) 

So the mean wall shear decays as 1 /V/ to an ultimate positive wall shear of mean 'AA2 with second harmonic. 
Thus the forcing in the Stokes layer is met by a particular solution confined to the Stokes layer but slipping at 
the wall. That generates vorticity / a homogeneous solution to prevent slip whose mean diffuses and spreads all 
the way out to the outer flow in time. 

The mean shear is just the time average of they integral of the forcing, whilst as eqn 5.13.19 [4] the streaming 
mean velocity is the time average of the y moment of the forcing. Now the Vj forcing is strong and near the wall 
whilst the Wi forcing is weak but further from the wall so its streaming mean velocity can be much more 
significant than its shear. 

The constant values in (11) at t=0 multiply homogeneous solutions, so these are exactly cancelled by their 
counterparts in wq h , and all that briefly survives are the /asymptotes due to these initial values. (This did not 
arise in the power of V/ basis as all / integrations are automatically zero at /=0.) 

The highest error term in w ip comes from the highest integrated term fo fo gdy of Q.\ multiplying d y g to 
give a term at the wall as cos(/- l An)Nnt whose homogeneous solution is unfortunately unknown. How the 
running mean of this boundary value decays should be how its mean diffused value would decay. The running 


www.ijmret.org 


ISSN: 2456-5628 


Page 44 








International Journal of Modern Research in Engineering and Technology (IJMRET) 
www.ijmret.org Volume 3 Issue 2 II February 2018. 

total is J f cos(t+ l An)Nt = (FresnelS-FrcsnclC)\2\t which tends to 0 so there is no mean transient 1 It response 
to add to those from the homogeneous pure harmonics (that do not start as cos’s). So then 

W\ v ~y€ sh:1 cos(t-yN2+V4K)s\m-2€ > ' l ' 1 sm(t-yN2)cost-V2€ ril sm(2t-yN2)+Vie rn cos(ylj2)-Vw ril sm{2t-^2y) 


W\h= -7e' v sin(2r-y)/4+ Vi erfc^ +13Q.|/8 


(26) 


Here the mean of trip’s d v g logdy’ is twice as strong at y=0 as the mean of Vi‘s half square V 2 UJ gdy] 2 . 
Thus each wall mean just comes from one in- phase “square” term , and wq p ’s is clearly if suprisingly twice as 
strong as Vi p ’s. Equation (13) is more direct than (26) for some of the shear terms. 


d v rv,=(5/2-7\2/4) sin(2t+'/47i) - l /2Q.,-13Q.,/8 =+.025sin(2f+‘/4jr) +W <tn -13/ 16V?tf 3 (27) 


The second harmonic component is about l/10 th of Vi’s , as expected from section I. Three terms in nq p , have 
ultimate mean shears but they cancel. Thus the ultimate mean w j flux into they=0 sink vanishes so there is no 
net mean nq p forcing integrated over y. The higher steady nq streaming comes from more of the mean uq 
produced further from the wall sink than Vi diffusing to infinity. The nq mean forcing from -dygJogdy} 
decays as e‘ r ' whereas the Vi mean forcing from -g 2 decays as e'" 2 so the former has twice they moment and 
so twice the steady streaming. 3 


d y (v! +w,) «. 353+.292sin(2f+ 1 /47t) +W V/ n -9/S^lnt 3 (28) 




Figure 4 Perturbation wall shears computed for G=sin t from t=0 by Taylor series. NB different shear scales. 

3 [4] twice attributes steady streaming under longitudinal acceleration to just the wq momentum transfer by mean 
udyW but calculates the net from \q+ ttq. Though a 50% enhanced nq is indeed solely responsible for the phase 
steady streaming in - 3 AWT{d x W+Wdxy) under a wave coreflow WRe{e ,t+ ' r }in Batchelor’s (5.13.20). 
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The combined shear is practically d y Vi except the transient shear from the outward diffusion of the steady 
streaming is trebled. In Taylor series calculations for G= sin t from t= 0 convergence was accelerated by the 
Pade method[6,7] . Fig 4 shows the high steady streaming transient shear of w { declining slowly as t~ Vl to the 
low ultimate second harmonic shear. The transient first harmonics in nq decline as t 2 . The large ultimate 
mean and second harmonic dominate the small steady streaming transient in d v V]. 

The support of a Commonwealth Scholarship is gratefully acknowledged. 

VI. Conclusions 

The perturbation of a time-diffusive boundary layer for longitudinal acceleration is the sum of a 
pressure gradient perturbation large near the wall and a ‘indraft’ perturbation with a modest peak in the middle 
of the boundary layer. The first gives an order of magnitude greater shear, but the second has twice the steady 
streaming flow which spreads outwards as the complementary error function of the diffusive similarity variable. 
This gives twice the transient shear declining as the inverse square root of time. The next order transient 
outward diffusion is an inverse time function of the diffusive similarity variable. For moderate time, a Taylor 
series in time numerically converges for almost 5 cycles of a harmonic. For any such (root) time series, the 
general particular solutions found are easily matched with exact homogenous ones to complete an explicit exact 
perturbation solution. 
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